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Abstract. Three dimensional continuous and discrete Fourier-like transforms, based on the three 
simple and four semisimplc compact Lie groups of rank 3, are presented. For each simple Lie group, 
there are three families of special functions (C-, S-, and i?-functions) on which the transforms are 
built. Pertinent properties of the functions are described in detail, such as their orthogonality within 
each family, when integrated over a finite region F of the 3-dimensional Euclidean space (continuous 
orthogonality), as well as when summed up over a lattice grid Fm C F (discrete orthogonality). The 
positive integer M sets up the density of the lattice containing Fm- The expansion of functions given 
either on F or on Fm is the paper's main focus. 
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1. Introduction 

New n-dimensional C- and S- and iS-transforms were recently described in [8-10]. Each transform 
is based on a compact semisimple Lie group of rank n and comes in three versions: analogs of Fourier 
series, Fourier integrals, and Fourier transforms on an n-dimensional lattice. They are named C- and S- 
and ^-transforms [19] in recognition of the fact that they can be understood as generalizations of the 
one dimensional cosine, sine, and exponential Fourier transform. 

The aim of this paper is to set the grounds for the 3-dimensional exploitation of transforms - described 
here as continuous transforms in a finite region F of the 3-dimensional Euclidean space R 3 , and also as 
discrete transforms of functions given on a lattice grid of points Fm C F of any density - in ready-to-use 
form. A positive integer M specifies the density. In some cases, the density of the grid is more flexible as 
it is dictated by not one, but two or three positive intergers. The grid could thus be made denser along 
certain axes. The symmetry of the lattice is dictated by the shape of F, or equivalently, by the choice of 
the Lie group. 

There are seven compact semisimple Lie groups of rank 3: 

SU(2) x SU(2) x SU(2) , SU(3) x SU{2) , 0(5) x SU{2) , G{2) x SU{2) , 

SU(3), 0(7), Sp(6). 

Throughout the paper we identify these cases by symbols that are often used for their respective Lie 
algebras: 

AxxAixAi, A 2 x Ax , C 2 x A ± , G 2 x A x , A 3 , B 3 , C 3 . 

The immediate motivation for this paper is our anticipation of the extensive use of the transforms 
given the need for processing the rapidly increasing amount of 3D digital data gathered today. In 2D, 
our group transforms offered only in some cases more than a marginal advantage, having emerged when 
satisfactory practical methods had already been developed and adequately implemented. So far, practical 
use of the functions in 2D rested on the fact that the continuous extension of the transformed lattice data 
displayed remarkably smooth interpolation between lattice points [3] (see also references therein). 

Special functions, which serve as the kernel of our transform (we call them C-, S-, and E'-functions or 
orbit functions), have simple symmetry property under the action of the corresponding affine Weyl group. 
The affine group contains as a subgroup the group of translations in R™, which underlies the common 
Fourier transform. This is the primary reason for the superior performance of our transforms, although 
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detailed comparisons, rather than examples, will have to provide quantitative content to substantiate 
such a claim. 

Other properties of the C-, S-, and ©-functions are not less important. 

Within each family, functions are described in a uniform way for semisimple Lie groups of any type and 
rank. In this work, we illustrate this uniformity by considering all seven rank 3 group cases in parallel. 
The price to pay for the uniformity of methods is having to work with non-orthogonal bases which are 
not normalized. 

The functions are defined in R™ and have continuous derivatives of all degrees. Their orthogonality, 
when integrated over the finite region F appropriate for each Lie group, was shown in [18]. The discrete 
orthogonality of C-functions in Fm has already been described in [16] and extensively used (see for 
example [4] and references therein) . The completeness of these systems of functions directly follows from 
the completeness of the system of exponential functions. 

A Laplace operator for each Lie group is given in a different set of coordinates. The C- and ^-functions 
are its eigenfunctions with known eigenvalues. On the boundary of F, the C-functions have a vanishing 
normal derivative, while S'-functions reach zero at the boundary. 

The functions have a number of other useful properties, which can be found in [8-10]. For example, 
the decomposition of their products into sums, the splitting of functions into as many mutually exclusive 
congruence classes as is the order of the center of the Lie group, etc. . . . 

A different but valid viewpoint on some of the special functions presented here, namely, functions 
symmetrized by the summation of constituent functions over a finite group [15], may turn out to be 
rather useful. The finite group, in the case of C- and .Si- functions, is the Weyl group of the corresponding 
semisimple Lie group. In the case of .©-functions, it is the even subgroup of the Weyl group. The Weyl 
group of SU (n) is isomorphic to the group S n of the permutation of n elements. This led to the recent 
implementations in [11,12], where instead of the Weyl group of SU(n), the S n group is used, and variables 
are given relative to an orthonormal system of coordinates. Furthermore, the even subgroup of S n is the 
alternating group. Related transforms were introduced most recently in [13, 14]. 

The paper is organized as follows. 

In Section 2, necessary definitions and properties of Lie groups and algebras are given and discussed. 
Semisimple Lie groups of rank 3 are considered in detail in Section 3. For each of these groups, we lay 
down the information necessary to construct and use their orbit functions for 3D continuous and discrete 
transforms. Section 4 is devoted to C-, S- and E- orbit functions and to their pertinent properties. 
Continuous and discrete orbit-function transforms are presented in Section 5. Some problems and possible 
applications arising in connection with orbit functions are formulated in the conclusion. An example of 
the application of orbit function transforms in the case of the group SU (2) x SU (2) x SU (2) is given at 
the end of the paper. 

2. Pertinent properties of Lie groups and Lie algebras 

The notion of orbit function of n variables depends essentially on the underlying semisimple Lie group 
of rank n. This section is intended to recall some of the standard properties of semisimple Lie groups/Lie 
algebras in general, and particularly those of rank 3, as well as properties of related Weyl groups. We also 
fix notation and terminology. Additional information about such Lie groups can be found for example 
in [2,5,7,20]. 

2.1. Definitions and notations. 

Let M™ be the real Euclidean space spanned by the simple roots of a simple Lie group G (equivalently, 
Lie algebra). The basis of the simple roots is hereafter referred to as the a-basis. An a-basis is not 
orthogonal and comprises simple roots of at most two different lengths. If a semisimple G is not simple, 
the a-bases of its simple constituents are pairwise orthogonal. 

For important practical (i.e. computational) reasons, it is advantageous to introduce also the basis 
of fundamental weights, hereafter referred to as the aj-basis. Moreover, for Lie groups with simple roots 
of two different lengths, it is useful to introduce bases dual to a- and w-bases, denoted here as a- and 
w-bases respectively. Occasionally it is also useful to work with the orthonormal basis {ei, e2, . . . e n } of 
K". Each subsection contains an explicit elaboration of these bases for the groups we consider. 
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The Cartan matrix C of G provides, in principle, all of the information needed about G. It is an 
n x n-matrix. In particular, it provides the relation between a-and w-bases: 

a = Civ lo = C~ 1 a , 

where 

' i,j€ {1,2,. ..,»}. 

Here (• , •) denotes the inner product in R". The length of the long simple roots is determined by an 
additional convention 

(c^long, Gtong) = 2. 

The dual bases are fixed by the relations 

Wi= K^)' <«*>«*> = = 

where <5^ is the Kronecker delta. 

The root lattice Q and the weight lattice P of G are formed by all integer linear combinations of the 
a-basis and w-basis, 

Q = Zai + Za 2 + • • • + Za n , P = 1mj\ + Zu 2 + ■ ■ ■ + Zui n . 

In general, Q C P, but in rank 3 Lie groups the equalities do not occur. Similarly, we can introduce the 
dual lattices Q and P. 

In the weight lattice P, we define the cone of dominant weights P + and its subset of strictly dominant 
weights P ++ 

P D P + = V + • • • + Z^"uj n D P ++ = Z> V + • • • + Z >0 uj n . 

For any simple Lie group G there is a unique highest root £ 

£ = rniai + m 2 Q!2 H h m„a„ = qidti + q 2 a 2 H h g„a„ (1) 

Coefficients mi and g,, i = l,n, are natural numbers, referred to as marks and comarks respectively 
They are well known for all simple Lie groups (see for example [5,7]). 

2.2. Weyl groups and their orbits. 

The Weyl group W(G) of a semisimple Lie group G is the finite group generated in R™ by reflections 
in n — 1 dimensional mirrors (hypcrplancs) orthogonal to the simple roots of G and containing the origin 
of R™. For a simple root «j, i = 1, n the corresponding reflection r ai is given by 

riX~r a .x = x—- — ■ — l —cti = x — (x, di)a», x £ R n . (2) 
{ati,ai} 

There is a general method for building Weyl group orbits. In Section 3, wc limit ourselves to recording 
the result of its application for all Lie groups of rank n = 3. 

It is assumed that we have fixed a semisimple Lie group of rank 3 and that we consider its weight 
lattice P. A VF-orbit can be generated from any point (a, b, c) <E R 3 , but in this paper we are almost 
always interested in W-orbits of points in P. It is convenient to specify the orbit by its unique point 
(a, b, c) = auj\ + buj 2 + CC03 with positive integer coordinates a > 0, b > 0, c > 0. We denote such a generic 
orbit by W(a, b, c). 

Every group contains the trivial one point orbit (0,0,0). The number of points of a generic orbit 
\W(a, b, c)| is equal to the order \W\ of the Weyl group. 

2.3. Affine Weyl groups and their fundamental domains. 

Consider the reflection with respect to the hyperplane containing the origin and orthogonal to the 
highest root £, see (1) 

We extend the set of n reflections r ai , given in (2), generating the Weyl group W, by one reflection r 

2£ 

tqx = r^x + £, where £ = , x e R™. 
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The resulting group transformations of R™, generated by n + 1 reflections ro, ri, . . . , r n , is referred to 
as the affine Weyl group W a ^ . The order of W a '' is infinite. 

The fundamental region F(G) C K" for any W aff (G) is the convex hull of the vertices {0, ^, . . . , ^}, 
where qi, i = l,n are comarks from (1), or equivalently, 

F(G) — {0, — , — , . . . , ^-}, where mu i = l,n are marks, see (1). (3) 
mi m 2 m n 

Repeated reflections of F(G) in its n — 1-dimensional sides results in tiling the entire space R™ by 
copies of F. We define grid Fm C F, depending on an arbitrary natural number M as given in (1), 

Fm= j^^i + H^ + '-' + ^w™ | «!,...,«„ Y^Simi<M > j . 

The number of points of the grid Fm is denoted by \Fm\, arid the volume of the fundamental region 
F with respect to the euclidian measure is denoted by \F\. 

Remark 1. In case G = G\ x G 2 , where G\ and G 2 are simple, the fundamental region of G is the 
Cartesian product of the fundamental regions of Gi and G 2 . The same holds for grids Fm, but for each 
simple constituent of G, the numbers Mi nd M 2 could be chosen independently. 

In the case of two different numbers Mi and M 2 , such that M\ mod M 2 = 0, the corresponding grids 
are related as follows Fm 2 C Fm x ■ 

2.4. Even subgroup W e of the Weyl group. 

Elements of the subgroup W e C W are formed by an even number of reflections which generate W. 
Each G— and S— function is built on a single Weyl group orbit of point A G P, and each E— function is 
built on an orbit of the even subgroup of the Weyl group, defined as follows. 

Consider a subgroup of W generated by an even number of reflections 



W e - {r n n 2 . . . r l2k \k G N, r H £ W, 1 = 1, 2k}. (4) 

W e is a normal subgroup of the index 2 of the Weyl group, i.e. 2\W e \ = \W\. Let us denote the orbit of 
point A G P with respect to the action of W e by We (A), and the size of this orbit by |W e (A)|, then 



m - / W e( X ) u W e (nX), for some r, G W, when A G P ++ 

( j " \ We(X), when A G P+ \ P++ 



(5) 



and one of the following relations holds true |W e (A)| = ||W(A)|, when A G P++ or |W e (A)| = |W(A)|, 
when A G P + \ P ++ . Note that each orbit of the even subgroup W e contains exactly one point from 

P e ■= P+ U nP++ , n e W . 

Similarly, we define the even affine Weyl subgroup 

Wf f = { ril r i2 . . . r l2k \k G N, r H G W aff , I = T^k}, 
and its fundamental region 

F e = FU r,F, (6) 

where F is a fundamental region of W a '' (3) and ti G ro, r\, . . . , r n . 

The same formula holds true for the grid on the fundamental region F e 

F eM = F M unF M , n ew aff . 

Remark 2. As it follows from (6), the fundamental region of W e is not unique, hence it can be chosen 
in such a way as to be convenient for a given application. 

Example 1. Consider the rank one compact simple Lie group SU{2) (the corresponding Lie algebra 
is A{). 

The Cartan matrix in this case is the lxl matrix G = (2). The root system consists of two roots ±a. 
The root and weight lattices are formed by integer multiples of the simple root and integer multiples of 
the fundamental weight u>. 

Q = {La}, P = {Lto}, where a — Cuo so that u) = \ot. 

Therefore P = QL)(Q + uj). 
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The Weyl group W of A\ is of order 2. It is a reflection group acting in R. We have W = {1, — 1} and 
W e = {1}. Consequently, a Weyl group orbit containing the point x ^ 0, also contains the point —x, but 
the orbit of the even subgroup of W consists of a single point, either x or —x. 

The fundamental region F is the segment with endpoints F = {0,w}. 

The grid Fm C F is fixed by the positive integer M and it consists of M + 1 points 

TP _ rn 1 2 M-i ii 

-TM I ' if ' M ' " ' ' M ' 1 J- 

3. Semisimple Lie algebras of rank three 

In this section, we provide specific information about compact semisimple Lie groups of rank 3, namely 
SU(2) x 5*7(2) x SU(2), 0(5) x 5?7(2), 5C7(3) x S77(2), G 2 X SU(2), and the simple groups Sp(6), 0(7), 
and 5f7 (4) . For each of the groups ,we lay down all of the information necessary to construct and use 
their orbit functions for 3D continuous and discrete transforms. 

Note, that only the weight lattices of the Lie groups £17(2) x SU(2) x SU(2), 0(5) x SU(2), 0(7) 
and Sp(6) display cubic symmetries. 

Below we present the Dynkin diagram for each semisimple Lie group of rank 3. On these diagrams, 
long and short simple roots a are respectively denoted by unfilled and filled circles and comarks are 
presented over the circles. 



3.1. The Lie algebra Ax X A t X A x . 




Figure 1. a) the fundamental region F of the Lie algebra Ax x Ax x Ax', x, y and z 
indicate respectively the orthogonal directions of the orthonormal basis {ei, 62, e^}; b) 
the grid F 2 ,x,3(Ax x Ax x A x ). 



This case is a straightforward concatenation of three copies of A\ (see Example 1). The Dynkin 
diagram and Cartan matrix with its inverse are the following, 

OOO c-(JSJ), c- = i(J!8). 

ax 0L 2 «3 

All simple roots of the same length equal to \[2. The bases of simple roots and fundamental weights 
are thus related by 

ai = 2uJi, Wj = |aj, &x = ai, &x = a>i, £=1,2,3. 
These bases are easily written in the orthonormal basis {ei, ei, &z\ 

014 = V2ei, u)i = -±=ei, Ui — ^e,, £=1,2,3. 
The highest roots for the simple subgroups A\ are given by the formulas 
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The fundamental region is a cube with vertices, see Fig. la 

F(Ai x A\ x Ai) = {0, uji, Ui + ujj, oji+ lu 2 + u 3 }, where j < i £ 1,2, 3. 
The volume of the fundamental region is given by the formula 

\F(A 1 x A x x A x )\ = • \oj 2 \ ■ |w 3 | = -^=- 

The grid Fm,m',m" C F is fixed by the independent choice of three positive integers M,M' and M", 
and consists of all the points 

F m ,m',m"( a i x M x Ar) = {^wi + 7^2 + ^3 I Sl < M, s[ < M', s'( < M'; si.ai.s? e Z^ } . 

The grid is cubic if M = M' = M", otherwise it is rectangular. The freedom to use unequal values of 
M, M' and M" may prove rather useful in the analysis of data given on rectangular, but not cubic grids. 
The number of points in the grid Fm,m',M" equals to 

\F m ,m',m»(Ai x Ai x Ai)| = • |F M <(.4i)| • |F M "(^i)| = (M + 1)(M' + 1)(M" + 1). 

Example 2. 

Consider the case M = M' = 1, M" = 3. There are 16 points of Fi ; i i3 (^i x A x x Ai). Explicitly, 
we have the following sets of integers [si, s' l7 s"] and the corresponding grid points in the w-basis (in this 
case the basis is orthogonal), (j^, jjp, jpr)- 

[0,0,0] -(0,0,0), [0,0,1] = (0,0,1), [0,0,2] = (0,0, |), [0,0,3] = (0,0,1), 

[0,1,0] = (0,1,0), [0,1,1] = (0,1, 1), [0,1,2] = (0,1, |), [0,1,3] = (0,1,1), 

[1,0,0] = (1,0,0), [1,0,1] = (1,0, i), [1,0,2] = (1,0, |), [1,0,3] = (1,0,1), 

[1,1,0] = (1,1,0), [1,1,1] = (1,1,1), [1,1,2] = (1,1, |), [1,1,3] = (1,1,1). 

Next, consider the case M = 2, M' = 1, M" = 3, see Fig. lb. The 16 points of F ljh3 {Ax x^x A x ) 
form a part of F 2t 1,3 (-A 1 x A\ x Ai). Note that they correspond to another set of [si, s' l7 s"]. More 
precisely, i*i,i,3 is the subset of ^2,1,3 with an even value of s\. Thus ^2,1,3(^1 x^x Ax) contains an 
additional eight points: 

[1,0,0] = (i, 0,0), [1,0,1] = (1,0,1), [1,0,2] = (1,0, §), [1,0,3] = (i, 0,1), 
[1,1,0] = (i, 1,0), [1,1,1] = (1,1,1), [1,1,2] = (1,1,|), [1,1,3] = (1,1,1). 

The Weyl group orbit of the generic point auj\ + bui 2 + CW3, a,b,c > 0, always consists of the eight 
points 

W {aM {A x xAiX Ax) = {(±o, ±6, ±c)}. 

Orbit sizes for arbitrary points are given by the relations 

\W( a ,b,c) I = 8, \W (aA0 ) I = |W (o ,o,c) I - \W (0M I = 4, 
1^(0,0,0)1 = 1. I WW) I = 1^(0,6,0)1 - \W i0 , , c) \ = 2. 



3.2. The Lie algebra A 2 x A x . 

The Dynkin diagram and Cartan matrix with its inverse are the following, 

cm^o Mi?!), C -'=»(!1I)- 

ol\ a 2 a 3 

Hence all the simple roots of the same length equal to \ 2. 

The bases of simple roots and fundamental weights are thus related by 

ot\ — 2ui - u 2 , u>i — |ai + 1q!2, 
a 2 = -wi + 2u 2 , u 2 = g«i + 3Q2, 
a 3 = 2ui 3 ; ui 3 = ^a 3 , 





= Oil, 




= Wl 


a 2 


= a 2 , 




= to 2 


a 3 


= «3, 


0)3 


= W3 
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Figure 2. a) the fundamental region F of the Lie algebra A 2 X A\\ x, y and z in- 
dicate respectively orthogonal directions of the orthonormal basis {ei, e 2 , 63} and the 
fundamental region of A 2 lies in the plane x + y + z = 0; b) the grid F2^(A 2 X Ai). 



In order to visualize the implied geometry, it is useful to represent the u> and a bases in the orthonormal 
basis. We have 

ax = (1, -1, 0) = e\ - e 2 , lox = (§, -\, -\) = fei - \e 2 - §e 3 = u>i, 

a 2 = (0, 1, -1) = e 2 - e 3 , lo 2 = (|, -|) = |ei + |e 2 - §e 3 = w 2) 

"3 = ^C 1 ' X ) = vf (ei + 62 + e3 )' W3 = 75 C 1 ' !) = 7f( ei + 62 + e3 ) = ^ 
The highest roots for the simple subgroups A 2 and Ax are given by the formulas 

( = ai+a 2 , £ = a 3 . 

The fundamental region F is a cylinder (see Fig. la) with an equilateral triangle as its base and w 3 
as its height. Its vertices are 

F(A 2 x Ai) = {0, lox, 0J2, U3, wi + oj 3 , lu 2 + w 3 }. 

The volume of the fundamental region is given by the formula 



\F{A 2 xAx)\ = \{pxM^) = \-\-l-^ 



2 -1 -1 
11-2 
111 



1 

2V6 



The grid Fm,m> {A 2 X Ax) C F(A 2 x Ax) is fixed by the independent choice of the two positive integers 
M and M'. It consists of all the points 

F m ,m-{A 2 x A x ) = [§w x + §& 2 + ^d> 3 I «i +s 2 < AT, s'j < M') sx,s 2 ,s' 1 e Z^ } . 

The number of points in the grid Fm,m'(A 2 x Ax) equals to 

\F m ,m>(A 2 x Ax)\ = \F M (A 2 )\ ■ \F M ,(Ax)\ = ((M + l) 2 - M*±ll) (M' + 1). 

Example 3. Consider the case M = 2, M' = 3, see Fig. 2b. There exist 24 points of F 2 , 3 

\F 2 ^A 2 xAx)\= ((2 + l) 2 - -2±ii) (3 + 1) = 24. 

Explicitly, we have the following sets of integers [sx , s 2 , ] and the corresponding grid points ( Q , |J , -^7 ) 
in cD-basis: 



[0,2,0] = 


(0,1,0), 


[0,2,3] 


= (0,1,1), 


[0,2,2] 


= (0,1, |), 


[0,2,1] 


= (o,i,|), 


[2,0,0] = 


(1,0,0), 


[2,0,3] 


= (1,0,1), 


[2,0,2] 


= (i,o,|), 


[2,0,1] 


= (i,o,|), 


[0,0,0] = 


(0,0,0), 
(0,1,0), 


[0,0,3] 


= (0,0,1), 

= (o,|,i), 


[0,0,2] 


= (0,0, |), 


[0,0,1] 


= (o,o,|), 


[0,1,0] = 


[0,1,3] 


[0,1,2] 


= (0,|,|), 


[0,1,1] 


= (0, 2 ,|), 


[1,0,0] = 


(5,0,0), 


[1,0,3] 


= (2,0,1), 


[1,0,2] 


= (- -) 

\ 2 ' ' 3>' 


[1,0,1] 


= (|,o,|), 


[1,1,0] = 


( 2 , 2' 0), 


[1,1,3] 


= (- - 1) 

\ 2 i 2 ' /' 


[1,1,2] 


= (---) 

V2' 2' 3/' 


[1,1,1] 


= (- - -) 

V2' 2' 3^ 



MARYNA NESTERENKO 1 ' 2 AND JlRl PATERA 



The Weyl group orbit of the generic point au)\ + bu 2 + cuj 3 , a, b, c > 0, consists of twelve points 
W {aM ((A 2 xii)) = {(a, 6, ±c), (-a, a + 6, ±c), (a + b, -b, ±c), (6, -(a + b), ±c), 

(-(a + 6), a, ±c), (-6, -a, ±c)}. 
Orbit sizes for arbitrary points are given by the relations 



\W (aAc) \ = 12, 
\W( a>0 ,o)\ = 3, 



l^(o,6,0)l = 6, 
1^(0,6,0)1 =3, 



\w. 



(o,0,c) 
(0,0, c) 



= 6, 
= 2, 



\w. 



(0,6,c) 
(0,0,0) 



= 6, 
= 1. 



,3. The Lie algebra C 2 x A 1 

y 





Figure 3. a) the fundamental region F of the Lie algebra C 2 x A\\ x, y and z indicate 
respectively orthogonal directions of the orthonormal basis {ei, e 2 , e 3 } and the funda- 
mental region of C 2 is the isosceles right-angled triangle in the plane z = 0, its right-angle 
is at the point i&i; b) the grid F 3 ^ 2 (C 2 x A%). 



The Dynkin diagram and Cartan matrix with its inverse are the following, 
2 1 1 



G 



2-10 
-2 2 
2 



1 Vo o 1/ 



mzo o 

oti a 2 a 3 

Hence o>i is the shorter of the simple roots. The relative lengths of the simple roots are set 

= 1 and (a 2 ,a 2 ) = (a 3 ,a 3 ) = 2. 
The bases of simple roots and fundamental weights are thus related by 

ai = 2cJi — uj 2 , oji — 1 J " ■ 

a 2 = —2uj\ + 2u 2l u> 2 
ot 3 = 2uj 3 ; uj 3 

In the orthonormal basis these bases have the form 



ai + a 2 , 

\oi 3 \ 



Oil 


= 2a X) 




= 2u) X 


&2 


= "2, 




= 0J2, 


&3 


= & 3 \ 


U) 3 


= UJ 3 . 



at = (0,1,0) = e 2 , 
a 2 = (1,-1,0) =ei- e 2 , 
a 3 = (0,0, y/2) = V2e 3 ; 



ui = (ii°) 



2 e l + 2 

W2 = (l,0,0) = ei, 
u 3 = (0,0,-4=) = 4=e 3 ; 



e2, 



LU 2 
U 3 



(1,1,1), 

(1,0,0), 
(0,0,^) 



\/2' V2 

The highest roots for the simple subgroups C 2 and Ai are given by the formulas 

£ = 2«i + a 2 , £ = a 3 . 
The fundamental region F(C 2 x 4i) is a cylinder (see Fig. 3a) with a triangular base. Its vertices 

F(C 2 x Ax) = {0, a> 2 , w 3 , |tl>i + w 3 , w 2 + ^3}- 

The volume of the fundamental region is given by the formula 

\F(C 2 xA 1 )\ = $([a2,& 1 ],S> a ) 



1 J_ 

4 ' s/2 



10 
110 
1 



1 

4V2 
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The grid Fm,m> C F is fixed by the independent choice of the two positive integers M and M'. It 
consists of all the points 

F m ,m'( c 2 x A t ) = {jfccbx + fa& 2 + ^7W 3 | 2 Sl + s 2 <M, si < M'; s u s 2 , s[ € Z^ } . 
The number of points in the grid i*M,M'(C 2 X Ai) equals to 

|F m ,m'(C 2 x A x )| = |F M (C 2 )| • |F M /(^)| = ([f ] + 1) (M + 1 - [f ]) (M' + 1), 
where [ • ] denotes the integer part of a number. 

Example 4. Consider the case M = 3, Af' = 2. There exist 18 points of ^3,2 (C2 x Ai), see Fig. 3b 

|F 3 , 2 (C 2 x A0| = ([§] + 1) (3 + 1 - [§]) (2 + 1) = 18. 

Explicitly, we have the following sets of integers [si,s 2 ,s' 1 ] and the corresponding grid points in the 
(i-basis, (jfc, ||, jjt): 



[1,0,0] = 


(|,o,o), 

(0,1,0), 


[1,0,1] = 


(|,o, i), 


[1,0,2] = 


(1,0,2), 


[0,1,0] = 


[0,1,1] = 


(o,ii), 


[0,1,2] = 


(0,1,2), 


[0,2,0] = 


(o,f,o), 


[0,2,1] = 


(0 2 1) 

v u ' 3 ' 2 /' 


[0,2,2] = 


(0,|,2), 


[1,1,0] = 


(3, 3, 0), 


[1,1,1] = 


(I I I) 

V3' 3' 2/' 


[1,1,2] = 


(3, 3, 2), 


[0,0,0] = 


(0,0,0), 


[0,0,1] = 


(0,0,1), 


[0,0,2] = 


(0,0,2), 


[0,3,0] = 


(0,1,0), 


[0,3,1] = 


(0,1,1), 


[0,3,2] = 


(0,1,2). 



The Weyl group orbit of the generic point au>x + buj 2 + cu>3, a,b,c > 0, consists of the following set of 
points 

W / (a : b,c)(C , 2 x Ax) = {±{a,b,c), ±(a,b, -c), ±(-o,a + 6,c), ±(a + 26, -6, c), ±(a + 26, -(a + 6), c), 
±(-a,a + b,-c), ±(a + 2b,-(a + b),-c), ±(a + 26, -6, -c)}. 

Orbit sizes for arbitrary points are given by the relations 

|W (o , 6 , c) | = 16, |W (0|6i o)| = 8, |W (a , . c) | = 8, 1^(0,6,0)1 = 8, 
\W (a , , )\ = i, |W(o,6,o)| = 4, |W(o,o,c)| = 2, |W(o,o,o)l = l- 



3.4. The Lie algebra G 2 X ii. 





Figure 4. a) the fundamental region F of the Lie algebra G 2 X A\\ x, y and z indicate 
respectively orthogonal directions of the orthonormal basis {ex, e 2 , e 3 } and the funda- 
mental region of G 2 is the half of the equilateral triangle in the plane z — 0; b) the grid 
F 4 , 3 (G 2 x Ax). 
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The Dynkin diagram and Cartan matrix with its inverse are the following, 

Oi\ a 2 Ol3 

Hence a 2 is the shorter of the simple roots. The relative lengths of the simple roots are set as 
{0-2,0-2) = § and (ai,ai) = (a 3 ,a 3 ) = 2. 

The bases of simple roots and fundamental weights are thus related by 

ai = 2u>i — 3lu 2 , u>i = 2ai + 3a 2 , o;i = 0\, uj\ = ui\, 
a 2 = —uJi + 2oj 2 , ui 2 = a 1 + 2a 2 , a 2 — 3a 2 , uj 2 = 3ui 2 , 
o 3 — 2w 3 ; w 3 = |a 3 ; a 3 = a 3 ; w 3 = a> 3 . 

Relative to the orthonormal basis, we have 

ai = (V2,0,0) = V2 ei , Wl = (-^,^,0) = ^ei + ^e 2) *i = (^, $,0), 

^ = (-75-75,0) = -^! + ^, lo 2 = (0,^,0) = ^e 2 , d&2 = (0,v^,0), 

a 3 = (0,0,V2) = V2e 3 ; c^ 3 = (0, 0, ^) = ^ e 3 ; w 3 = (0,0,^). 

The highest roots for the simple subgroups are given by the formulas 

£ = 2ai+3a 2 , £ = a 3 . 
The fundamental region F is a cylinder (see Fig. 4a) with a triangular base. Its vertices are 

F(G 2 x Ai) = {0, |w 2) w 3 , ±uj 1 +lj 3 , |w 2 +w 3 }. 

The volume of the fundamental region is given by the formula 



\F(G 2 x A 1 )\ = I^J.d&s) = I • 575 • 75 • 7^ 



1 v^3 

V20 

1 



where [•, •] denotes the vector product. 



The grid Fm,m> C F is fixed by the independent choice of two positive integers M and M'. It consists 
of all the points 

Fmm>{G 2 x Ai) = {^wi + |Jw 2 + ^uo 3 I 2si + 3s 2 < M, s[ < M'; s 1 ,s 2 ,s' 1 e Z^ } . 
The number of points in the grid Fm,m' equals to 

\F M ,M"(G2 x - |F M (G 2 )| • 1^(^)1 - ^[f ] + 1 + E j ( M ' + X )> 

where [ • ] is the integer part of a number. 

Example 5. Consider the case M — 4, M' — 3. There are 16 points of F^ 3 (see Fig. 4b) 

( [l] \ 

|F 4 , 3 (G 2 x A 1 )\ = [§] + 1 + J2 J (3 + 1) = (1 + 1 + (2 + 0)) • 4 = 16. 

Explicitly, we have the following sets of integers [si,s 2} s' 1 ] and the corresponding grid points in the 
w-basis, (&,&,;&): 

[0,0,0] = (0,0,0), [0,0,1] = (0,0, |), [0,0,2] -(0,0, |), [0,0,3] = (0,0,1), 

[0,1,0] = (0,±,0), [0,1,1] = (0,±,±), [0,1,2] = (0,1,1), [0,1,3] = (0,±,1), 

[1,0,0] = (|,0,0), [1,0,1] = (±,0,±), [1,0,2] = (i, 0,|), [1,0,3] = (±, 0,1), 

[2,0,0] = (§, 0,0), [2,0,1] = (|,0,i), [2,0,2] = (1,0, |), [2,0,3] = (±, 0,1). 

Next, consider the case M = 3, M' = 2. F 3 , 2 (G 2 x Ai) consists of the following points: 

[0,0,0] = (0,0,0), [0,0,1] = (0,0, \), [0,0,2] = (0,0,1), 
[1,0,0] = (i, 0,0), [1,0,1] = (i,0,i), [1,0,2] = (i 0,1), 
[0,1,0] = (0,|,0), [0,1,1] = (0,±,|), [0,1,2] = (0,1,1). 
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In this case, none of the nine points of F 3i2 (G2 x Ax) coincides with a point of ^4,3 (G?2 x Ax). This is due 
to the fact that lattice densities M — 3 and M' — 2 do not correspondingly divide the densities M — 4 
and M' = 3. 

The Weyl group orbit of the generic point au)\ + bu>2 + cuj 3 , a, b, c > 0, consists of 24 points 
W( a>btC )(G 2 x Ai) = {±(a,b,c), ±(-a,3a + 6,c), ±(o + 6,-6,c), ±(2o + 6, -(3a + 6),c), 
±(-(a + 6), 3a + 26, c), ±(-(2a + 6), 3a + 26, c), 

±(a, 6, — c), ±(— a, 3a + b, — c), ±(a + b, — b, — c), ±(2a + 6, — (3a + 6), — c), 
±(-(a + 6), 3a + 2b, -c), ±(-(2a + 6), 3a + 2b, -c)}. 
Orbit sizes for arbitrary points are given by the relations 

\W {aAc) \ = 2A, \W (am \ = 12, \W (afiiC) \ = 12, |W (0>M | = 12, 
1^,0,0)1 = 6, 1^(0,6,0)1 = 6, |W (0) o,c)l = 2, |W ( o,o,o)l = l. 




Figure 5. a) the fundamental region F of the Lie algebra A3; x, y and z indicate respec- 
tively orthogonal directions of the orthonormal basis {ex, e 2 , e%\, and the fundamental 
region of A 3 is the tetrahedron; b) the grid F 3 (A 3 ). 



The Dynkin diagram and Cartan matrix with its inverse are the following, 
Oi\ a 2 a 3 

Hence, all simple roots of the same length equal to \[2. 

The bases of simple roots and fundamental weights are thus related by 

ai = 2a>i— W2, wi = |ai + |a 2 + |a 3 , a, = on, i£ {1,2,3}. 

a 2 = -oj\ + 2uj 2 - £j 3 , w 2 = |ai + a 2 + 2" a 3, u>i = u)i, i€ {1,2,3}. 
a 3 = -W2 + 2w 3 ; w 3 = I ax + \a 2 + §a 3 ; 

Relative to the orthonormal basis, we have 

a x = (1,-1,0) =e x - e 2 , 
"2 = (0, 1, -1) = e 2 - e 3 , 
a 3 = (3, 3, |) = + e 2 + 4e 3 ); 
The highest root £ is given by the formula 

£ = ctx + a 2 + a 3 . 



^2 = (I, I, -|) = |(2ei + 2e 2 - e 3 ) = w 2 , 
w 3 = (5, |, 5) = 2"( e i + e 2 + e 3) = &3- 
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The fundamental region is a pyramid (see Fig. 5a) with vertices 

F(A 3 ) = {0, Cbi, uj 2 , w 3 }. 
The volume of the fundamental region is given by the formula 

\F{M)\ = WiM^) = \-\-k-\ 

where [•, •] denotes the vector product. 

The grid Fm C F is fixed by the choice of one positive integer M. It consists of all the points 

F M (A 3 ) = + ffW2 + ffW3 I Si +S2 + S3 < M] «i,«2,S3 € Z^°} . 

The number of points in the grid Fm equals to 

1 M 

\F M (A 3 )\ = -J2(M + l-i)(M + 2-i). 

i=0 

Example 6. Consider the case M = 3. There are 20 points of F 3 (A 3 ) (see Fig. 5b). Explicitly, we have 
the following sets of integers [si, s 2 , S3] and the corresponding grid points in the w-basis ( jfe, , ||): 



[0,0,0] 


= (0,0,0), 


[0,0,3] = 


(0,0,1), 


[0,3,0] 


= (0,1,0), 


[3,0,0] = 


(1,0,0), 


[2,0,1] 


= (- i) 

V3' ' 3/' 


[2,1,0] = 


(1 1 0) 

V3, 3, u jj 


[1,0,2] 


= fl- 2) 

l3' u ' 3^' 
= (1.0,0), 


[0,1,2] = 


(o,| 


1) 


[1,2,0] 


= (3 ' 3 ' 0), 


[0,2,1] = 


(0, 3, 3), 


[1,0,0] 


[0,0,1] = 


(0,0, 


1), 


[0,1,0] 


= (0,|,0), 

= a i i) 

V3' 3' 3'' 


[1,1,0] = 


(3,3,0), 


[1,0,1] 


= fi i) 

I3, u , 3 h 


[0,1,1] = 


(o,| 
(0,1 


1) 


[1,1,1] 


[2,0,0] = 


(1,0,0), 


[0,0,2] 


= (o,o,|), 


[0,2,0] = 


0) 



The Weyl group orbit of the generic point au\ + bui 2 + c^ 3 , a,b,c> 0, consists of 24 points 
W( aAc )(A 3 ) = {(a, 6, c),(-a, a + b, c),(a + b, -b, b + c),(a, b + c, -c),(b, -(a + b), a + b + c), 
(-a, a + b + c, -c),(-(a + b), a, b + c), (a + b, c, -(b + c)), 
(a + b + c, -(b + c), c),(b, c, -(a + b + c)), (b + c, -(a + b + c), a + b), 
(-b, -a, a + b + c),(-(a + b), a + b + c, -(b + c)), (-(a + b + c), a, b), 
(a + b + c, -c, -b),(b + c, -c, -(a + b)),(-b, b + c - (a + b + c)), 
(c, -(a + b + c), a), (-(b + c), -a, (a + b)),(-(a + b + c), a + b, -b), 
(c, -(b + c), -a), (-(b + c), b, -(a + b)),(-c, -(a + b), a),(-c, -b, -a)}. 
Orbit sizes for arbitrary points are given by the relations 

\W {aAc) \=24, \W {aA0) \ = 12, \W {a . a . c) \ = 12, 1^(0,^)1 = 12, 
1^,0,0)1=4, |W(o,6,o)l=6, |W( ,o, c )|=4, |W(o,o,o)| = l. 

3.6. The Lie algebra B 3 . 

The Dynkin diagram and Cartan matrix with its inverse are the following, 

^-i(?il)- 

ax a 2 a 3 

Hence ot\ is the shorter of the simple roots. The relative lengths of the simple roots are set as 
(ai,ai) = 1 and (a 2 ,a 2 ) = (a 3 ,a 3 ) = 2. 

The bases of simple roots and fundamental weights are thus related by 

«i = 2cji — 0J2, lui = ai + a 2 + a 3 , a\ = a\, oj\ = wi, 

a 2 = —ijj\ + 2lu 2 — 2u 3 , lu 2 = ct\ + 2a 2 + 2a 3 , a 2 = a 2 , uj 2 = lu 2 , 
a 3 = -lu 2 + 2uj 3 ; lo 3 = \a\ + a 2 + \a 3 ; a 3 = 2a 3 ; lu 3 = 2lo 3 . 

Relative to the orthonormal basis, we have 

ai = (l,-l,0) = ei-e 2 , wi = (1,0,0) = e u u>i = (1,0,0) = e u 

a 2 = (0,1,-1) = e 2 -e 3 , uj 2 = (1, 1, 0) = e 1 + e 2 , ui 2 = (1, 1, 0) = ei + e 2 , 

a 3 = (0,0, 1) = e 3 ; lj 3 = (5, \, \) = \(a + e 2 + e 3 ); w 3 = (1, 1, 1) = e x + e 2 + e 3 . 
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Figure 6. a) the fundamental region F of the Lie algebra B 3 ; x, y and z indicate 
respectively orthogonal directions of the orthonormal basis {ei, 62, 63} and the face 
defined by vertices {0,wi, hoJ 2 } lies in the plane z — 0; b) the grid Fi(B 3 ). 



The highest root £ is given by the formula 

£ = oti + 2a 2 + 2a 3 . 
The fundamental region is a pyramid (see Fig. la) with vertices 

F(B 3 ) = {0, (&!, iw 2 , §& 3 }. 
The volume of the fundamental region is given by the formula 



\F(B i )\ = \{[<hMM = \-\-\ 



100 
110 
111 



= ■57, where [•,•] denotes the vector product. 



The grid Fm C F is fixed by the choice of an integer M. It consists of all the points 



F M (B 3 ) = { 



M 



M 



U) 2 



f-3 



Si 



2s 2 + 2s 3 <M; sus 2 ,S3 G Z- } 



The number of points in the grid Fm equals to 

\FM(B 3 )\ = ([f]+l)([f] [M±l]+M + 1 



^[f]) 



i=0 



where [ • ] is the integer part of a number. 

Example 7. Consider the case M = 4, see Fig. 6b. There are 14 points of F4 

[|] 

\F 4 (B 3 )\ = ( [f] + 1) ( [f ] [*±1] + 4 + 1 - i±2 [§] ) + ^ = 3(2 • 2 + 4 + 1 - 3 • 2) + + 1 + 4 = 14. 

i=0 

Explicitly, we have the following sets of integers [si, 52,53] and the corresponding grid points in the 
w-basis, (jfe, ||, H): 



[0,0,0] -(0,0,0), 
[0,2,0] = (0,|,0), 
[0,0,4] = (0,0,1), 
[0,1,0] = (0,±,0), 



[0,0,1] = (0,0, \), [0,0,2] 

[o,i,i] = (o,i,i), [0,1,2] 

[0,0,3] = (0,0, |), [1,1,0] 
[1,0,0] = (i, 0,0), 



(o,o,|), 
(o,l,|), 
(iio), 



[1,0,1] = (±,0,±), 
[1,0,2] = (±,0,±), 
[2,0,0] = (i, 0,0). 



The Weyl group orbit of the generic point au)\ + bu 2 + cuj 3 , a, 6, c > 0, consists of 48 points 
W(o,6,c)(£3) = {±(0, 6, c),±(-a, a + 6, c),±(o + 6, -6, 26 + c),±(a, 6 + c, -c), 

±(6, -(a + 6), 2a + 26 + c),±(-a, a + & + c, -c), ±(-(a + b), a, 26 + c), 
±(a + &, 6 + c, -(26 + c)),±(a + 6 + c, -(6 + c), 26 + c), 
±(6, a + b + c, -(2a + 26 + c)),±(6 + c, -(a + b + c), 2a + 26 + c), 
±(-6, -a, 2a + 26 + c),±(-(a + 6), a + 26 + c, -(26 + c)), 
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±(a + 2b + c, -(b + c), c),±(-(a + b + c), a, 2b + c), 
±(a + 6 + c, b, -(26 + c)),±(a + 26 + c, -(a + 6 + c), c), 
±(6 + c, a + 6, -(2a + 26 + c)),±(-6, a + 26 + c, -(2a + 26 + c)), 
±(b + c, -(a + 26 + c), 2a + 26 + c),±(-(a + 26 + c), a + 6, c), 
±(a + 26 + c, -6, -c),±(-(6 + c), -a, 2a + 26 + c), 
±(-(a + 6 + c), a + 26 + c, -(26 + c))}. 
Orbit sizes for arbitrary points are given by the relations 

|WW)I=48, |WW,o)l=24, |W (o>0lC) | = 24, |W (0Ac) | = 24, 
\W (am \ = 6, |W (0lM )| = 12, |W (0 ,o,c)l = 8, |W(o,o,o)l = 1- 



3.7. The Lie algebra C 3 . 




a 



Figure 7. a) the fundamental region F of the Lie algebra C 3 ; x, y and z indicate 
respectively orthogonal directions of the orthonormal basis {ei, e%, e 3 } and the face 
defined by vertices {0, \Co 2 } lies in the plane z = 0; b) the grid F±{C 3 ). 



The Dynkin diagram and Cartan matrix with its inverse are the following, 

w w — w V -2 2 / 2 \ 2 4 3, 

«i a 2 a 3 

Hence a 3 is the longer of the simple roots. The relative lengths of the simple roots are set as 
(ai, ai) = (a 2 , a 2 ) = 1 and (a 3 , a 3 ) = 2. 

The bases of simple roots and fundamental weights are thus related by 

a\ = 2loi — u 2 , Wi = eti + a 2 + |a 3 , «i = 2ai, U\ = 2a>i, 

a 2 = — cjj + 2t^2 — w 3 , W2 — a i + 2a2 + 0:3, a 2 = 2a 2 , a; 2 = 2co 2 , 
a 3 = — 2lo\ + 2lu 3 ; uj 3 = a\ + 2a 2 + |a 3 ; a 3 = a 3 ; ii 3 =aj 3 . 

In the orthonormal basis these bases have the form 

= (72 ' -73 • °) = 71^ - e ^ ^ = ( A' °' 0) = 7f ei ' = °' 0) ' 

a2 ^(°'72^7f) = 7f( e2 - e 3)' W2 = (7^72' ) = 7f( e i + e2 )' ^ 2 = (V2,V2,0), 

a 3 = (0, 0, V2) = V2e 3 ; W3 = (^, ^, ^) = ^( ei + e 2 + e 3 ); ^ 3 = (^, ^, ^). 

The highest root £ is given by the formula 

£ = 2«i + 2a 2 + a 3 . 

The fundamental region is a pyramid (see Fig. la) with vertices 

F(C 3 ) = {0, |d>2, ci 3 }. 
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The volume of the fundamental region is given by the formula 

\F(C 3 )\ = \([\u; u \u 2 ],u 3 ) = \-\-\- V2-V2-± 

where [•, •] is the vector product. 

The grid Fm C F is fixed by the choice of the integer M. It consists of all the points 

F M (C 3 ) = {^loi + j^Lu 2 + j^Lb 3 | 2ai + 2s 2 + s 3 < M; s 1; s 2 , s 3 e Z^ } . 

The number of points in the grid Fm equals to 

m 

i^m(c 3 )i = ([f ] + 1) ([f ] m +m+ i - m±i m 

i=0 

where [ • ] is the integer part of a number. 

Example 8. Consider the case M = 4. There are 14 points of F^(C 3 ), see Fig. 7b. 

[I] 

i^(c 3 )i = ([§] + 1) ([|] m + 4 + 1 - m + e * 2 = 14. 

i=0 

Explicitly, we have the following sets of integers [si,s 2 ,S3] and the corresponding grid points in the 
w-basis, (f?' 17' m ) : 

[0,0,0] - (0,0,0), [0,0,1] = (0,0, 1), [0,0,2] = (0,0, 1), [0,0,3] = (0,0, 1), 

[0,1,0] = (0,|,0), [0,l,l] = (0,i,i), [0,0,4] = (0,0,1), [2,0,0] = (§, 0,0), 

[1,0,0] = (i 0,0), [1,0,1] = (1,0,1), [1,0, 2] = (1,0, |), [1,1, 2] = (1,1, |). 

[1,1,0] = (1,1,0), [1,1,1] = (1,1, |), 

The Weyl group orbit of the generic point au>i + bu 2 + cu> 3 , a,b,c > 0, consists of 48 points 
W {aAc) (C 3 ) = {±(a, b, c), ±(-a, a + b, c), ±(a + b, -b, b + c), ±(a, 6 + 2c, -c), 

±(6, -(a + 6), a + 6 + c), ±(-a, a + b + 2c, -c), ±(-(a + 6), a, b + c), 

±(a + b, b + 2c, -(b + c)), ±(a + b + 2c, -(6 + 2c), 6 + c), 

±(6, a + fc + 2c, -(a + b + c)), ±(b + 2c, -(a + b + 2c), a + b + c), 

±(-b, -a, a + b + c), ±(-(a + b), a + 2b + 2c, -(b + c)), 

±(a + 2b + 2c, -(6 + 2c), c), ±(-(a + 6 + 2c), a, b + c), 

±(a + b + 2c, b, -(b + c)), ±(a + 2b + 2c, -(a + b + 2c), c), 

±(b + 2c, a + b, -(a + b + c)), ±(-b, a + 2b + 2c, -(a + b + c)), 

±(b + 2c, -(a + 2b+2c), a + b + c), ±(-(a + 2b + 2c), a + b, c), 

±(a + 2b + 2c, -b, -c), ±(-(b + 2c), -a, a + b + c), 

±(-(a + b + 2c), a + 26 + 2c, -(b + c))}. 

Orbit sizes for arbitrary points are given by the relations 

|W (0iM |=48, |W (Oi6i 0)|=24, |W (0i o, c) |=24, |W (0 , 6 , c) | = 24, 
1^(0,0,0)1=6, |W(o,6,o)| = 12, |W( ,o,c)| = 8, |W(o,o,o)| = l- 

4. Orbit functions 

In this section, we define what we mean by C-, S- and .E-functions, specified by a given point A e Z™. 
We also show some of the properties inherent to those functions. Namely, the following properties 
are of interest, their pairwise orthogonality (using the appropriate scalar product, an integral over the 
fundamental region), their discrete orthogonality (again, using a properly defined scalar product, a sum 
over the discrete grid), their product can be represented as a sum, and they are eigenfunctions of the 
Laplace operator. 
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4.1. Definitions, symmetries and general properties. 

We start with the C-functions. The C-function C\{x), A e P + is defined as 

C x (x):= e 2wi ^' x \ ieR", (7) 
new* 

where W\ is the Weyl group orbit generated from A. 

If in (7) we restrict ourselves to the orbit of the even subgroup W e \, then we define E- function E\{x), 

AeP e 

E x {x) := e 2 * iM , (8) 

The definition of an S- function 5a (x), A e P ++ is almost identical, but the sign of each summand is 
determined by the number of reflections p(^t) necessary to obtain \i from A 

S x (x) := (-l) pM e 2 ™ M (9) 

Of course the same n can be obtained by different successions of reflections, but all routes from A to 

\i will have a length of the same parity, and thus the salient detail given by , in the context of an 
S-function, is meaningful and unchanging. 

For different families of orbit functions, the A (represented in the cj-basis) are taken from different 
sets, namely 

A e {Z^ V + V + Z^°w 3 }, for C -functions; 

A e {Z^ V + Z^ 2 + Z^°lu 3 } U ri{Z> V + Z> V + Z >0 w 3 }, neW for E -functions; 

A e {Z> V + Z >0 w 2 + Z >0 w 3 }, for S -functions. 

In particular, this implies that, for S-functions, the number of summands always equals to the size of 
the Weyl group. 

In the case of x e Fm (the coordinates of x are rational), the C-, S- and iJ-functions are formed by 
roots of unity. Therefore, there can only be a finite number of possible orbit functions that can take 
distinct values on the points of Fm, these functions are orthogonal on the grid. The number of pairwise 
orthogonal orbit functions on Fm coincides with the size of Fm, including the boundary in the case of 
C- and ^-functions and excluding the boundary in the case of , S-functions. 

Note that in the 1-dimensional case, C-, S- and _E-functions are respectively a cosine, a sine and an 
exponential functions up to the constant. 

All three families of orbit functions are based on semisimple Lie groups of finite order, the number of 
variables coincides with the rank of the corresponding Lie algebra. 

In general, C-, S- and S-functions are the finite sums of exponential functions, therefore they are 
continuous and have continuous derivatives of all orders in R™. 

It is easy to prove that C- and S'-functions are invariant with respect to the action of both Weyl W 
and affine Weyl W a ^ groups (see e.g. [8,9]) and that ^-functions are invariant with respect to the action 
of W e and W° ^ groups (see e.g. [6]). Therefore it is enough to consider them only on the fundamental 
domain of their affine Weyl symmetry groups. 

The S'-functions are antisymmetric with respect to n — 1-dimensional boundary of F . Hence they are 
zero on the boundary of F. The C-functions are symmetric with respect to n — 1-dimensional boundary 
of F. Their normal derivative at the boundary is equal to zero (because the normal derivative of a 
C-function is an S- function) . A number of other properties of orbit functions are presented in [8-10] . 

4.2. Calculation of scalar products. 

Here we present the rules necessary for the calculation of the scalar products of the vectors given in 
the different bases. 

Let vectors u = [u\, U2, ■ ■ ■ , u n ) = u\uji + U2LO2 + ■ ■ ■ + u n uj n and v = viu>i + 1*2^2 + ■ ■ ■ + v n uj n are 
represented in the w-basis, then their scalar product is calculated as the matrix product 

n n , . 

(u,v) = uCv 1 = ^ UiVj{uJi,tjjj) = ^ u i v j — 3 ^ 3 
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here Ci , = ^f 1 ^ C~} and C is the Cartan matrix. 

If vectors x = (xi, x 2 , ■ ■ ■ , x n ) = X\a\ + x 2 a 2 + • • • + x n a n and y = y\a\ + y 2 a 2 + ■ ■ ■ + y n a n are 
represented in the a-basis, then their scalar product is calculated as the matrix product 

n n I \ 

(x,y) = xCy 1 = x tVj { ail a 3 ) = ^ .*-,//, J ' ' <",,„. 

here Ci,j — dj and C is the Cartan matrix. 

Square lengths of the simple roots {aj,aj) are always indicated in the Dynkin diagrams of the semisim- 
ple Lie groups, see e.g. [5,7]. 

It is also useful to remind that 

(ai 7 cbj) = (&i,u)j) = (ei,ej) = Sij, where Sij is the Kronecker delta. 

4.3. C-, S-, and _E-functions as eigenfunctions of the Laplace operator. Consider the functions 
C\(x), E\(x) and S\{x) and suppose that the continuous variable x is given relative to the orthogonal 
basis. In the case of Lie algebra A n we use orthogonal coordinates xi,x 2 , . . . ,x n +i and coordinates 
x\,x 2 , ■ ■ ■ ,x n for B n , C n and D n (the orthogonal bases for these algebras are well known and can be 
found e.g. in [8]). 

The Laplace operator in orthogonal coordinates has the form 

d 2 d 2 d 2 

A = — h 7; !-■■■ + 7; — , where k = n (or k = n + 1 for A n ). 

ox i ox 2 oxk 

For the algebras A n , B n , C n and D n , the Laplace operator gives the same eigenvalues on every 
exponential function summand of an orbit function with eigenvalue — 47r(A, A). 

Hence, the functions C\{x), E\{x) and S\{x) are eigenfunctions of the Laplace operator: 

C\(x) \ ( C x {x) 

A ( E x (x) ^-4^ 2 (A,A) E x (x) 
Sx(x) J V S *( x ) 

Now we consider the continuous variable x given relative to the w-basis. Let A denote the Laplace 
operator, where the differentiation d Xi is made with respect to the direction given by Wj. 

" C 

A = - — l -^—rd Xi d Xi , where C is the Cartan matrix. 

^ {ati,oti} 

1,3=1 

It is known in Lie theory that the matrix of scalar products of the simple roots is positive definite, 
moreover our definition makes matrix ^ a symmetric, hence it can be diagonalized and the Laplace 
operator could be transformed to the sum of second derivatives by an appropriate change of variables. 

Thereby, using the results of Section 3 we can write the explicit forms of the Laplace operators given 
in the aj-basis for all semisimple Lie algebras of rank 3 

(d 2 xl +d 2 X2 +d 2 3 , for A 1 xA 1 xA 1 ; 

d 2 xl -d Xl d X2 +d 2 X2 +d 2 3 , for A 2 xA i; 

2d 2 Xl -2d Xl d X2 +d 2 X2 +d 2 X3 , for C 2 xA i; 

3d 2 Xi -5d xl d X2 +d 2 2 +d 2 X3 , for G 2 xA 1 ; 

d Xl - d xl d X2 + d 2 X2 - d X2 d X3 + d 2 3 , for A 3 ; 

9 Xl - d Xl d X2 + dl 2 - 2d X2 d X3 + 2d 2 3 , for B 3 ; 

[2d 2 Xl ~2d Xl d X2 +2d 2 X2 -2d X2 d X3 +d 2 X3 , for C 3 . 

4.4. Continuous orthogonality. 

For any two squared integrable functions <j>(x) and ip(x) defined on the fundamental region F, we 
define a continuous scalar product 

(<Kx),il>(x)) ~ J p <t>(xj$(x)te. (10) 

Here, integration is carried out with respect to the Euclidean measure, the bar means complex conjugation 
and x € F, where F is the fundamental region of either W or W e . 



A = <^ 
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Any pair of orbit functions from the same family is orthogonal on the corresponding fundamental 
region with respect to the introduced scalar product (10), namely 

(C\(x),Cy(x)) = \W X \ ■ \F\-5m, (11) 

(S x (x),S x ,(x)) = \W\-\F\-5 xx ,, (12) 
(E x (x),E x ,(x)) = \W eX \-\F e \.S xx ,, (13) 

where 5 XX > is the Kronecker delta, \W\ is the size of Weyl group, \W X \ and \W e x \ are the sizes of Weyl 
group orbits, and \F\ and |-F e | are volumes of fundamental regions. All of the necessary information for 
each semisimple Lie algebra of rank 3 can be found in Section 3. In particular, \F e \ = 2\F\ and \W e x \ is 
defined by formula (5). 

Proof of the relations (11,12,13) follows from the orthogonality of the usual exponential functions and 
from the fact that a given weight /i€P belongs to precisely one orbit function. 

Therefore each family of orbit functions forms an orthogonal basis in the Hilbert space of squared 
integrable functions C? (F) . Hence functions given on F can be expanded in terms of linear combinations 
of C-, S- or i?-functions. 

4.5. Discrete orthogonality. 

Let us denote the discrete grid of the fundamental region as Fm in the general case, even though in 
some cases it is determined by more than one positive integer 

for C-functions, 
7 m = { Fm \ dF for S'-functions, 
for _E-functions. 

A discrete scalar product of two functions <p(x) and ip(x) given on Fm (including C-, 5- and _E-functions) 
is dependent on this grid and defined by the bilinear form 

\Fm\ 

((j}(x),ip(x)) M = ^2 e(x i )(j)(x i )ip(xi), Xi € F M - (14) 

»=1 

Here e(xi) is the number of points conjugate to Xi on the maximal torus of the Lie group. The value of 
e(xi) is given by the formula 

\W Xi \ for C-functions, 
s{xi) — { \W\ for S'-functions, 

| W e Xi | for E- functions. 

Again, as in the continuous case, the C-, S- and i?-functions are pairwise orthogonal, i.e. 

(C x (x),C y (x)) M = Yl \W x JCx(xijC^faj = \Wx\-\A M \-5xx,, (15) 

i=T^N 

(S x (x),Sx>(x)) M - \W\ ]T S x (x i )S x 4^) = \W\ ■ \A M \-6 XX >, (16) 
i=i7/v 

(Ex(x),E y (x)) M = \WexMx{xi)hAxi) = \W e x\-\A M \-5xx'- (17) 

i=IJV 

Here \W Xi \ and \ W e Xi \ denote the sizes of the orbits of the Weyl group and its even subgroup. 

Am is the IF- invariant Abelian subgroup of the maximal torus T of the simple compact group corre- 
sponding to W 

\Fm\ 

A m := {wx\x e F M , w G W} and \A M \ = ^ 1^1- 

»=i 

Proof of the orthogonality relations can be found in [18]. 

All necessary data for the computation of the coefficients and discrete scalar products is given in 
Section 3. 
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5. C-, S- AND S-TRANSFORMS 

For different fixed m G R n the set of exponential functions S^ e 2 ^ l ( m ^ x ) ; x S i?™} determines contin- 
uous and discrete Fourier transforms on W 1 . In much the same way, the orbit functions (which are a 
symmetrized version of exponential functions defined in Section 4) determine an analogue of the Fourier 
transform. 

In this section, we introduce the essentials of the continuous and discrete C-, S- and ^-transforms. 
The discrete transform can be used for the continuous interpolation of values of a function f(x) between 
its given values on a grid Fm- 



5.1. Continuous transforms. 

Each continuous function on the fundamental region with continuous derivatives can be expanded as 
the sum of C-, S- or _E-functions. Let f(x) be a function defined on F (or F e for E- functions) , then it 
may be written that 

f(x) = cxC x (x), c x = \W x \- 1 \F\- 1 (f(x),C x (x)); (18) 

AeP+ 

/(*)= £ c » S ^ c x = \W\- 1 \F\- 1 (f(x),S x (x)}; (19) 

AGP++ 

f{x) = ]T c x E x (x), c x = | We x \- 1 \F e \~ 1 (f(x),Ex(x)). (20) 

AGP e 

Here (•, •) denotes the continuous scalar product of (10). Direct and inverse C-, S- and ^-transforms of 
the function f(x) are in (18), (19) and (20) respectively. 



5.2. Discrete transforms. 

Let Am € P be the maximal set of points, such that for any two A, A' G Am the condition of discrete 
orthogonality holds for any of the families of orbit functions in ((15), (16) or (17)). 
Then we have the following discrete transforms for the function f(x): 

f(x) = bxC x (x), x e F M , b x = ; (21) 



AeA M 

f{x) = £ b x S x (x), xeF M , b x = ^' JVM ; (22) 



AeA a 



fix)= ^ b x E x ix), xeF eM , b x = //'^\ M ■ (23) 



AGAm 

Here (•, -)m denotes the discrete scalar product given of (14). 



(Cx,Cx)m 

if, S x )m 
(S x ,S x )m 

if, E X ) M 
(E x , E x )m 



5.3. Continuous extensions. 

Once the coefficients b x of the expansions (21) (22), (23) are calculated, discrete variables Xi in Fm 
may be replaced by continuous variables x in F 

fcont(x) := 52 5 aCa(2i), xeF; 
AeA M 

fcont(x) := 52 b ^S X (x), x e F; 

AeA M 

fcont(x) := 52 b * E *( X )' X G F <=- 

The function f con t(x) smoothly interpolates the values of f(x{), i = 1, 2, . . . , \Fm\- At the points x i7 we 
have the equality f con t(xi) = f(xi). 
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Remark 3. If we calculate more than one discrete transform on the same grid Fm and use the same set 
of orbit functions <p\(x), A S Am, then it is reasonable to pre-compute and save the matrix 

/ 0AiOl) I ) 

/ TX" X" — T - ■ ■ ■ TT" T. 



B 



4>x n (xi) 

(4>\ n ,4>x„)m 



</>a„Qw) 



7 



This would save valuable computation time, especially for large M, since the coefficients b\ are easily 
calculated as a matrix product b = B ■ / * . 

Moreover, the matrix B does not depend on the function that is to be expanded into series, therefore 
it need be calculated only once for each M and can be repeatedly used. 

6. Concluding remarks 

• Each of the transforms described here is based on a compact semisimple Lie group of rank 3. 
All seven types of such Lie groups were considered. Our goal was to provide the tools for the 
expansion of functions of 3 variables given on a bounded region D of an Euclidean space R 3 . 
The variables can be either continuous or discrete (lattice grid points). The symmetry of the 
lattice is the Weyl group of the Lie group. The bounded region D has to be scaled to fit into the 
fundamental region F of the corresponding Lie group. In case of functions given on a lattice grid, 
the scaling has to be accompanied with the matching density of the grid points in F. Fortunately, 
the formalism admits choosing any density one may need. The scaling resulting in the inclusion 
D C F is not unique. Various options may be considered for specific functions. 



























1 \ 



















Figure 8. One of the possible applications of orbit functions is the construction of 
unknown transitional data (smooth interpolation). Fig. 8 shows how additional frames 
could be added to a film. In much the same way, the continuous deformation of the 
picture can be proceed and a 3-dimensional image can be created from corresponding 
2-dimensional layers or cuts. 



• The uncommon special functions of our transforms are defined for compact semisimple Lie groups 
of any type and rank [19]. Their continuous and discrete orthogonality in F is assured [18]. Unlike 
the translation symmetry required in traditional Fourier expansions, the symmetry group of C-, 
S-, and _E-functions is the appropriate affine Weyl group, which contains the translations in M 3 
as a subgroup. 

• The uniformity of our approach, as to the type of the rank 3 Lie group, is illustrated here 
by considering the seven cases in parallel. The price paid for uniformity is the exploitation of 
non-orthogonal bases, a, u-bases and their duals whenever necessary. The majority of practically 
useful digital data usually given on cubic/square lattices with the simplest symmetry group. Only 
more costly experimental installations may use denser lattice arrangements of data collectors. 

It should be pointed out that, at least for one type of transform, it is possible to avoid paying 
the price i.e. of having to work with orthonormal bases. The Weyl group of SU(n) is isomorphic 
to the permutation group S n of n elements. Recently introduced transforms [11-14], based on 
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S n and on its alternating subgroup, exploit orthonormal bases in R n , although even there, the 
corresponding fundamental regions do not have orthogonal adjacent faces, in general. 

• There is an additional freedom of choice whenever the underlying Lie group is not simple. Suppose 
that group is a product of two simple Lie groups, G 1 x G 2 . The fundamental region is then the 
Cartesian product of ^(G 1 ) and F(G 2 ). 

For the expansion of class functions on the product group, we can combine G-functions on one 
with ^-functions on the other. Similarly, we can combine G- or ^-functions with iJ-functions, 
enlarging appropriately the fundamental region of the .E-functions. 

In much the same way, discretization can proceed differently on F{G 1 ) then on F(G 2 ). The 
corresponding integers M\, Mi that fix it can be as different as one desires. Thus the density of 
grid points in F(G 1 ) x F(G 2 ) may be very different on the two orthogonal components. 

• A number of other properties of the orbit functions may prove to be useful (see [8-10] and 
references therein) . Let us point out that each of the three types of functions split into mutually 
exclusive congruence classes. For a given semisimple Lie group, the number of congruence classes 
equals the order of the center of the Lie group. 

• The possibility to introduce the G-, S-, and E- functions by summation over a finite noncrystallo- 
graphic Coxeter groups instead of the Weyl group of a Lie group appears to be rather interesting. 
In 3D there is just one such group H 3 , the icosahedral group of order 120. Most of the proper- 
tics carry over to this case in a simple straightforward way. The exception is the orthogonality, 
continuous or discrete. There is an analog of the fundamental region, but no lattice. Its role, 
perhaps, should be played by some quasicrystal? 

Appendix A. Example of G-transform on A x x A\ x A\ 

As an example, we chose to interpolate a discretization of a known function, namely the Gaussian 
function shown in equation (24) 

g(x) = e-^-P) 2 = e -(*-p)-(*-p), x e r3 7 (24) 

where p e F(Ai x Ai x Ax) is a fixed point inside the fundamental region. 

The first test to be undertaken is to sample the function in the points of the grid Fm for several values 
of M. The continuous extension of g(x), calculated from points of Fm, is 

T(x) = 5>G A (x), with 6 A = ^1^L. (25) 
ftp {Cx ' Cx) 
For different values of M, Fig. 9a represents the error, as defined by the integral 

\T(x)-g(x)\dF. (26) 



IF 

The integral (24) was calculated using a simple Monte Carlo method with 10000 randomly chosen points. 
We made sure that the granularity of the randomly generated points was much higher than that of the 
grid F M . 

Fig. 9b shows the standard deviation of the sample set. It decreases in much the same way the error 
did in function of M . 

As a means to show what takes place visually, we give a plot of T as compared to g on the parametric 
line given by the transformation 

L:R^m. 3 , k^(k,k,k). (27) 

Fig. 10a and Fig. 10b show the comparison between T (L(k)) and g (L(k)) for k e [0, 1], with M = 4 and 
M = 10, respectively. As seen when comparing Figs. 10a and 10b, the degree of oscillation increases as 
M increases, but the total error decreases. This is exactly as one would expect. 
Consider two-dimensional cuts on the parametric surface given by 

SiK^IR 3 , (k,l) A (k + l,k + k,k + l) with k e [0, \] and Z e [0, §]. (28) 

A two dimensional cut from the graphs of g is presented in Fig. 11. Figs. 12a and 12b are two 
dimensional cuts from the graphs of T with M = 4 and M = 10, respectively. 
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Figure 9. a) the error (26) of the C-transform (25) as function of M; b) the standard 
deviation of the difference between the Gaussian function (24) and its interpolation (25) 
as function of M. 




FIGURE 10. Comparison of T (L{k)) (25) with g (L(k)) for k G [0, 1] (27), with M = 4 
(a) and M = 10 (b). 




Figure 11. Two dimensional cut (28) of the Gaussian function g (S(k, I)). 



In the case of Al x Al x Al (as well as in the other cases, albeit less trivially), the set of pairwise 
orthogonal functions is not unique, i.e. one can scale all A G P by integer multiplier s ■ (a, b, c) with 
(a, b, c) G Z 3 , s £ Z. This leads to a shifted system of orbit functions, that take the exact same values on 
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Figure 12. Two dimensional cut (28) of the interpolation T (S(k, I)) for M = 4 (a) and 
M = 10 (b). 

the discrete grid Fm, but do in fact differ in their continuous behavior. Such a shifted system of functions 
can be understood as higher harmonics of the original functions on the grid Fm- 

Examples with original and shifted systems of orbit functions follow. Fig. 13a shows g and T with 
M = 6, for T computed using the fundamental set of C-functions and Fig. 13b shows the same graph, 
but for a T that was computed using the set of C-functions shifted by a factor of (M, 0,0). 




a b 



Figure 13. Comparison between g(L(k)) and T (L(k)) (25) computed using the fun- 
damental set of C-functions (a) or computed using the shifted by a factor (M, 0, 0) set 
of C-functions (b) for k € [0, 1] (27), with M = 6. 

The result is enlightening. The case using the shifted set of orbit functions exhibits a degree of error far 
above what could be tolerated, due to the fact that the higher harmonics of the C-functions, as expected, 
oscillate more. This derogatory result should not, however, dismiss the use of such higher harmonics. 
The Gaussian function extrapolated here is very smooth, and should thus be extrapolated with a sum of 
smoother functions. But in real life applications, the data could be quite chaotic, thus the use of higher 
harmonics could be useful. 

It is also possible not to shift the entire set, but only a specific subset, and at that, not all should be 
shifted by the same factor. Careful consideration must be taken, because the functions do not obey the 
simplistic rule that a shifted function is equal to its counterpart. In fact, the set of adjacent functions is 
paired according to a reflective symmetry, but not according to translation symmetry. 
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